Quasiparticle Formation and Optical Sum Rule Violation in Cuprate Superconductors 
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Using a simple model for the frequency dependent scattering rate, we evaluate the in-plane optical integral 
for cuprate superconductors in the normal and superconducting states. In the overdoped region, this integral 
is conserved. In the optimal and underdoped region, though, the optical integrals differ, implying a lowering 
of the in-plane kinetic energy in the superconducting state. This sum rule violation, due to the difference of 
the non Fermi liquid normal state and the superconducting Fermi liquid state, has a magnitude comparable 
to recent experimental results. 
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In superconductors, there is a dramatic change in the 
conductivity due to opening of an excitation gap in the 
finite frequency response, and the formation of a zero 
frequency 6 function peak representing the dissipationless 
response of the condensate. The change is such as to 
preserve the optical sum rule, in that the finite frequency 
weight removed by the opening of the excitation gap is 
recovered by the condensate peak [Q| . 

In cuprate superconductors, though, there is experi- 
mental evidence that the sum rule is violated for c-axis 
conductivity Over the measured frequency range, 
more weight is present in the condensate peak than can 
I be accounted for by the loss of finite frequency weight. 
Since the total optical sum rule must be preserved, the 
extra weight in the condensate peak is coming from out- 
side this frequency range. This is unusual, since in clas- 
sic superconductors, the change in the optical integral 
is exhausted over a frequency range of order 4A, where 
A is the superconducting gap. Anderson has stressed 
that such sum rule violations are connected to the lack of 
quasiparticle poles in the normal state, and their emer- 
gence in the superconducting state. This is unlike the 
' premise of BCS theory, where quasiparticles are assumed 
to exist in the normal state. 

Although the sum rule violation for the c-axis response 
is profound, its contribution to the condensation energy 
is small due to the smallness of the c-axis kinetic energy 
in the cuprates. If kinetic energy effects are to play a 
role in the condensation energy, then they must be com- 
ing from the in-plane response, since the in-plane kinetic 
I energy is quite large, of order an eV Q. This largeness, 
though, means that the violation is difficult to see. That 
is, a 1 meV change in the kinetic energy would represent 
^1% change in the optical integral. Recently, though, 
two groups have claimed to observe such a change. El- 
lipsometry data on optimal and underdoped Bi2212 ||] 
have been quantified as corresponding to a change in the 
kinetic energy of 1 meV. The same kinetic energy change 
has been inferred from reflectance data |6| on an under- 
doped Bi2212 film, though no such change could be re- 



solved in an overdoped film. These results are intrigu- 
ing, since a 1 meV kinetic energy savings per plane is in 
excess of the condensation energy inferred from specific 
heat data 

In this paper, using a simple model for the frequency 
dependent scattering rate based on angle resolved pho- 
toemission (ARPES) and infrared data, we calculate the 
change in the optical integral from the normal to the su- 
perconducting state, and find its sign and magnitude to 
be comparable to these recent findings. 

The full optical integral, integrating over all energy 
bands, is proportional to the bare carrier density over 
the bare electron mass, and thus must be conserved. Of 
greater interest here is the optical response of the band 
around the Fermi energy, correlating with the experimen- 
tal data which are typically integrated out to an energy 
of order the plasma frequency (1 eV). This leads to a 
consideration of the single band sum rule H 
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where the restriction of a to the single band response is 
implicit, and where a is the in-plane lattice constant, V 
the unit cell volume, and 
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with N the number of k vectors, ek the bare dispersion as 
defined by the effective single band Hamiltonian |P, and 
rik the momentum distribution function. For a Hamil- 
tonian with near neighbor hopping jj], Ek is equivalent 
to minus the kinetic energy {Ekin = X]fc ^kn-k), but in 
general these two quantities differ. 

For free electrons, the inverse mass tensor is a constant 
in momentum, and thus this integral is conserved due to 
charge conservation. This is not generally the case, since 
the sum of the inverse mass tensor over the Brillouin 
zone vanishes ]lO| . When considering the change in this 
integral between different electronic states, the emphasis 
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in the past has been on a possible change in the inverse 
mass tensor In general, though, we expect to 

be invariant, and therefore the change should instead be 
due to changes in n^. A simple case is BCS theory [|, 
where the kinetic energy increases in the superconducting 
state due to particle-hole mixing. If a near neighbor tight 
binding model applied, the BCS optical integral would be 
smaller in the superconducting state than in the normal 
state, opposite to the recent experimental results. 

The BCS model, though, assumes the existence of 
quasiparticle poles in the normal state. It is straightfor- 
ward to demonstrate that the kinetic energy can indeed 
be lowered in the superconducting state if the normal 
state is a non Fermi liquid and the superconducting state 
a Fermi liquid ^ . This occurs if the effect of quasiparticle 
formation on sharpening Uk is larger than the smearing 
due to particle-hole mixing. This effect is anisotropic in 
momentum, due to anisotropies in the scattering rate and 
the d-wave order parameter. Given these anisotropies 
and the anisotropy of the mass tensor, it is not obvious 
what the effect of the kinetic energy lowering will be on 
the optical integral, since a near neighbor tight binding 
model is inadequate to describe ek- In addition, ARPES 
measurements indicate a substantial doping dependence 
of the scattering rate, which implies that the sum rule 
violation will also be doping dependent. 

We start by considering a simple model for the fre- 
quency dependent scattering rate, based on fits to 
ARPES data at the {n,0) point M]. This was used in 
work on the condensation energy |p| and the c-axis sum 
rule [ p^ . The model assumes a large frequency indepen- 
dent scattering rate in the normal state, consistent with 
the broad Lorentzian lineshapes. In the superconducting 
state, the broad peak is replaced by a sharp peak at the 
superconducting gap energy, followed at higher binding 
energy by a spectral dip, then a broad maximum (the 
"hump"). This change is modeled by cutting off ImS at 
the energy of the spectral dip. The resulting E is 
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where wq is the spectral dip energy. This self-energy is 
then used in the spectral function (l4| 
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where Z = 1 — S/cj. For this form of S, the spectral 
function has two S functions located at ±E, where E sat- 
isfies the pole condition (denominator of Eq. 4 vanishes). 
Such poles always exist for E < ojq because of the log 
divergence of ReE at ±ojo- The weight of the poles are 
determined as \dA~^{±E)/du!\. In addition, there 
are incoherent pieces for |a;| > uq. 

For now, we assume F is fc independent, ujq is also 
assumed to be k independent, as implied by ARPES ex- 
periments pM. efc is taken from a six parameter tight 



binding fit to normal state ARPES data |T^. For the 
order parameter, the d-wave form cos(fca;a) — cos(fcya) is 
assumed. The k sum is done using a 100 by 100 grid 
in the irreducible quadrant of the zone. The quasiparti- 
cle pole weight contribution to Uk is analytic The 
incoherent contribution is evaluated by trapezoidal in- 
tegration. We consider the T=0 limit, and thus Uk — 
j'^^A{uj)dLj. In practice, the lower cut-off is taken to 
be -10 eV. In the normal state with no lower cut-off, 
nk = l/2- tan^i (e/F) /tt. 
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FIG. 1. — Ai?K versus (a) isotropic scattering rate and (b) 
antinodal scattering rate (for various nodal scattering rates). 
tiJo=71meV and Am,ai=32meV. 

In Fig. la, IS.Ek = E^ — E^ (where N denotes normal 
state and S superconducting state) is plotted as a func- 
tion of F. The strong dependence on F is expected, since 
as F increases, the change in nt (Arifc = — nf ) be- 
comes increasingly pronounced, leading to a larger sum 
rule violation. This implies that the sum rule violation 
becomes larger as the doping decreases, since F from 
ARPES measurements increases with underdoping. 

One issue with Fig. la is the rather large value of the 
sum rule violation for realistic values of F (the antin- 
odal scattering rate from ARPES is ^100 meV for op- 
timal doping). It is known, though, that the scatter- 
ing rate from ARPES is a strong function of momentum 
p9[ . We consider a simple model for the anisotropy with 
Ffe — Tn[^ + CR{cos{kxa) — cos(fcya))2/4], where Fat is 
the nodal scattering rate and FAr(l -I- cr) the antinodal 
one. In Fig. lb, we plot AEk versus cr for various Fat, 
and find that it rapidly saturates with cr, and thus with 
the antinodal scattering rate. We have also considered 
the influence of the anisotropic pseudogap on the normal 
state ||2^, and found this had little effect on AEk- 

To gain further insight, we plot in Fig. 2a the in- 
tegrand of AEk as a function of momentum. Note 
that the overall integral is negative, with negative re- 
gions corresponding to unoccupied states near the d-wave 
node ((0,0) — (tTjTt) Fermi crossing) and occupied states 
near the (7r/2, 0) points, and positive regions to occu- 
pied states near the node and unoccupied states near 
the antinode ((tt, 0) — (7r,7r) Fermi crossing). To under- 
stand this, we plot two curves on Fig. 2a, one the Fermi 
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surface, the other the zero of the inverse mass tensor. 
In our model is equal to 1/2 on the Fermi surface, 
and thus An^ changes sign there. Therefore, the optical 
integrand, which is the product of the inverse mass ten- 
sor times Ari/t, changes sign each time one of these two 
curves is crossed. From this, one can easily understand 
the various sign regions in the plot. 



For LOc, fits to ARPES are consistent with a value of 500 
meV I2I. 




FIG. 2. (a) /S.Ek versus k (red is positive, green near zero, 
blue negative). The curves are the Fermi surface and zero 
of the inverse mass tensor, (b) A(Vnfe ■ Vet) versus k (red 
is positive, blue near zero). r=150meV, (j;o=71meV, and 
Amai=32meV. 

Perhaps more instructive is to convert the optical in- 
tegral to the equivalent one involving — Vn^ • Ve^ (using 
Greens theorem for periodic functions ||ic|]). The result- 
ing integrand is plotted in Fig. 2b, and as expected, is 
localized about the Fermi surface. The important point 
is that the integrand peaks at the node. This can be eas- 
ily understood. In the superconducting state there are 
quasiparticle poles, but at the node, A/j = 0, so there is 
a true step discontinuity in n/^ there. As one moves away 
from the node, Aj. increases from zero, and so iVn^l 
decreases in magnitude. From Fig. 2b, it is easy to ap- 
preciate the result of Fig. 1 that the optical integral is 
sensitive to the nodal scattering rate and not so sensitive 
to the antinodal one. As our model was motivated by 
fitting ARPES data in the antinodal region of the zone, 
this indicates that a model based directly on the nodal 
region should be considered. 

Normal state ARPES data are consistent with a 
scattering rate of the form -ImE = Ffe + ck|a;|, where F^ 
has the anisotropy described above and a is momentum 
independent |^l| . For simplicity, we will assume that 
both of these terms have an infrared cut-off at loq as 
we did for the F model. When determining ReS, it is 
important to provide an ultraviolet cut-off to ImS. A 
hard cut-off at ujc leads to a log singularity in ReS at loc- 
Rather, we take ImS to saturate at uJc- This gives 
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where S — Eq, + Sr. The normal state S is obtained by 
setting luq — 0. luq is the energy of the dispersion kink 
along the zone diagonal, which is the same energy as the 
spectral dip at (7r,0) |lq|, and thus loq is k independent. 
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FIG. 3. (a) -AEk and (b) Ek versus F. a=0.75, 
a'o=71meV, Amax=32meV , and ajc=500meV. 

For illustrative purposes, we consider first the case with 
no anisotropy in F. In Fig. 3, we plot the sum rule vi- 
olation versus F for a typical value of a. The variation 
with F is similar to Fig. la despite the presence of a 
substantial a term. Thus for a "pure" marginal Fermi 
liquid (F = 0), the sum rule violation is essentially zero. 
The reason is that with F = 0, the normal state posseses 
quite sharp spectral peaks, and thus the change in Uk 
when going into the superconducting state is reduced. 

To make quantitative comparisons to experiment, re- 
alistic values of F^ and a as a function of doping are 
needed. We can obtain them from the optics data. In 
Fig. 4a, we plot 1/t{oj) for four Bi2212 samples in the su- 
perconducting state extracted from reflectivity data . 
The linear high frequency behavior is of the form a + buj. 
Let us relate these parameters to F^ and a. The a term 
is easy to obtain, since it is k independent. At T=0, 1/r 
is an average of -2ImE over a frequency range of to a; 
|p3| . Since the a term is linear in uj, then a = b. 

The Ffc term is a different story. If it were isotropic, 
then F = a/2. For the anisotropic case, these two quan- 
tities are related by a Fermi surface integral. We can 
do this analytically by replacing the anisotropy term 
by cos^(2</)), where 4> is the Fermi surface angle (the 
node is at (/) = 7r/4). We find that the Fermi ve- 
locity along the Fermi surface can also be fit to the 
same anisotropic form. The resulting transport inte- 
gral is |15) 2r(0) = [/ #w(0)/F(</))]/ / d0?;((/i) where 
v{(j)) — vn[1 + viiCos'^{2(j))] {v is the modulus of the ve- 
locity) and F((/i) = FAr[l CRC0S^{2(j))]. From ARPES 
p9[ cr = 3, and from the tight binding fit, vr — —0.72. 
Solving, we find that F^^ = 3.4r(0). 

For the other parameters, we note that the deviation 
of 1/r (w) from linearity sets in at an energy A^ax + i^o, 
where ujq is the single particle scattering rate gap. This 
is easily shown from the Kubo bubble by dressing one 
of the two lines. From ARPES and tunneling, ujq = 
Amax + ^res , whcrc ujres IS the energy separation of the 
peak and the dip ||12,eJ,E5|. This is found to vary with 
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doping as 5Tc [g5| , and we use this to extract Amax from 
the optics scattering rate gap. The resulting Amax val- 
ues are consistent with ARPES Q and tunneling p^ ] 
measurements. 
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FIG. 4. (a ) 1/t{uj) versus to for various Bi2212 samples 
from Ref. |E3] (OD overdoped, OPT optimal doped, UD un- 
derdoped). (b) Calculated sum rule violation {—AEk) versus 
doping, X. The curve is Tc. The parameters (meV) extracted 
from (a) are Tn (1, 22, 27, 37), a (.65, .75, .88, .98), luq (54, 
71, 76, 83), and Amax (24, 32, 41, 54) for OD70, OPT90, 
UD82, and UD67 respectively. Also shown in (b) are the ex- 
perimental results (open squares from Ref. open diamonds 
from Ref. |5|). The theoretical doping trend in (b) is due to 
the increasing offset in 1/r seen in (a). 

These values are used to determine the sum rule viola- 
tion versus doping, shown in Fig. 4b. We find no sum rule 
violation for the overdoped sample. This is the expected 
BCS like behavior, and is consistent with the experimen- 
tal result on an overdoped film ||]. For the other samples, 
we find a sum rule violation which increases from 1.5 to 
2.2 meV as the doping decreases. The doping trend is 
consistent with the reported experimental results (also 
shown in Fig. 4b), although the values are perhaps too 
large by a factor of two. This may be due to the approx- 
imation of using a hard infrared cut-off on the scattering 
rate in the superconducting state. Still, given the sim- 
plicity of our model, and the substantial experimental er- 
ror bars, the agreement with experiment is surprisingly 
good. The doping trend in our model is due to the in- 
crease in r^v with undcrdoping. We also note that the 
kinetic energy change is about twice — AEk- 

As for where the extra condensate weight is coming 
from, we note that the the experimental optical integrals 
balance at an energy cut-off of about 2 eV This 
value is comparable to the Mott gap of the insulator, 
so we speculate that the extra weight comes from the 
upper and lower Hubbard bands. This would be in accord 
with the more delocalized nature of the electrons in the 
superconducting state. 

In conclusion, using a simple model for the frequency 
dependent scattering rate, we can understand recently 
reported results for the sum rule violation for the in- 
plane conductivity. The effect is due to the formation of 
quasiparticles in the superconducting state, and confirms 
earlier speculations by Anderson 0]. As the doping in- 



creases into the overdoped region, we find the sum rule 
violation goes away, consistent with the more Fermi liq- 
uid like nature of the normal state. 
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